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We calculate the network fluid regime and phase diagrams of binary mixtures of patchy colloids, 
using Wertheim's first order perturbation theory and a generalization of Flory-Stockmayer's theory 
of polymerization. The colloids are modelled as hard spheres with the same diameter and surface 
patches of the same type, A. The only difference between species is the number of their patches -or 
functionality-, f^ * 1 and f^ (with f^ > f^). We have found that the difference in functionality 
is the key factor controlling the behaviour of the mixture in the network (percolated) fluid regime. 
In particular, when /jj 2 > 2/^ the entropy of bonding drives the phase separation of two network 
fluids which is absent in other mixtures. This changes drastically the critical properties of the 
system and drives a change in the topology of the phase diagram (from type I to type V) when 
> 2. The difference in functionality also determines the miscibility at high (osmotic) pressures. 
If — fjp — 1 the mixture is completely miscible at high pressures, while closed miscibility gaps 

at pressures above the highest critical pressure of the pure fluids are present if — > 1. We 
argue that this phase behaviour is driven by a competition between the entropy of mixing and the 
entropy of bonding, as the latter dominates in the network fluid regime. 



I. INTRODUCTION 

Recent advances in the engineering of well-defined col- 
loidal particles, with designed surface patterning in the 
nanometer-to-the-micrometer range, open up the possi- 
bility of tailoring their behaviour at the macroscopic level 
[lHH. The resulting anisotropic or patchy particles have 
become the center of very active research, one focus be- 
ing the development of large-scale fabrication techniques 
which are required for the exploitation of the novel ma- 
terials in a range of applications Q. 

The current materials revolution based on patchy col- 
loidal particles results from recent breakthroughs in par- 
ticle synthesis and from the fact that these anisotropic 
colloids may be viewed as the molecules of future materi- 
als that can be designed to (self)assemble into functional 
structures. Indeed, the analogy between patchy colloids 
and molecules provides a powerful theoretical framework 
that exploits the anisotropy of the interactions to control 
the self-assembly of macroscopic structures Q. 
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Modeling and computational efforts describing patchy 
particle interactions and their macroscopic properties 
have followed the technological lead and produced a num- 
ber of interesting results The primitive model of 
patchy colloids consists of hard-spheres with / patches 
on their surfaces. Patchy particles attract each other if 
and only if two of their patches overlap. The attraction 
between particles is short ranged and anisotropic: The 
patches act as bonding sites and promote the appear- 
ance of well defined clusters, whose structure and size 
distribution depend on the properties of the patches (/ 
and the energy of attraction) and on the thermodynamic 
conditions (density and temperature). 

In the last 5 years Sciortino and co-workers inves- 
tigated, using Monte Carlo simulations, the phase be- 
haviour and the connectivity of the fluid phases of the 
primitive model of patchy colloidal particles, finding a 
large number of interesting properties. Among others, 
they established that /, the number of patches or bond- 
ing sites per particle, is the key parameter controlling 
the location of the liquid- vapour critical point (Tol . ITlj . 
They showed that, for low values of / (approaching 2), 
the phase separation region is drastically reduced and 
low densities and temperatures can be reached with- 
out encountering the phase boundary. These low den- 
sity ("empty") phases were shown to be network (per- 
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colated) liquids (see Fig. [TJ, suggesting that, on cool- 
ing, patchy particles with low functionality assemble into 
glassy states of arbitrary low density (gels). It was also 
shown that, at the critical point, it is always two network 
fluids that become identical. Very recently, Ruzicka and 
co-workers have reported the first experimental evidence 
of empty liquids in dilute suspensions of Laponite (l2l |. 
Their properties were found to be similar to those pre- 
dicted by the primitive patchy colloidal models. 

Remarkably, the results of the simulations of patchy 
colloidal particles are very well described by classical liq- 
uid state theories: Wertheim's first order perturbation 
theory fl3| predicts correctly the equilibrium thermody- 
namic properties; Flory-Stockmayer [l4| theories of poly- 
merization describe quantitatively the size distributions 
of the clusters of patchy particles, including the appear- 
ance of network (percolated) fluids. 

Here we extend the investigation of the thermodynamic 
properties and connectivity of patchy colloidal fluids to 
binary mixtures of patchy colloidal particles. As for pure 
fluids, we will use Wertheim's thermodynamic perturba- 
tion theory (in its extension to mixtures) [l5| and Flory- 
Stockmayer 's theory of percolation (generalized to mix- 
tures - see section HTB|) . 

In order to quantify the effect of decreasing function- 
ality (approaching 2) on the critical point of pure fluids, 
Sciortino and co-workers considered binary mixtures of 
particles with 2 and 3 patches However, they (delib- 
erately) ignored the entropy of mixing in the theoretical 
calculations. Instead, they used Wertheim's theory for 
pure fluids with a non-integer functionality set to the (av- 
erage) functionality of the mixture, assumed to be fixed 
by the functionality of the mixture at the critical point 

EE!- 

A few other studies reported results for the phase di- 
agram of polymer solutions, based on computer simula- 
tions and Wertheim-like theories (l6l . [l7| . Full critical 
lines connecting the critical point of the pure polymer 
to the critical point of the pure solvent were reported, 
and good agreement between the simulation and the- 
oretical results for the coexisting liquid densities and 
compositions was found at low pressures These 
models, however, are not directly related to the bi- 
nary mixtures of patchy particles referred above, as they 
consider monodisperse (fixed length) polymer chains in 
monomeric isotropic solvents. 

In what follows, we address, explicitly, the interplay 
between the entropy of mixing and the entropy of bond- 
ing in the network fluid regime of models of binary mix- 
tures of patchy particles. We restrict our study to mix- 
tures of patchy particles that differ by the number of 
bonding sites or patches on each species. 

We have found that the phenomenology in the network 
fluid regime is rich and somewhat surprising as the inter- 
play between the entropy of mixing and the entropy of 
bonding plays a dominant role, which is absent in other 
mixtures and/or regimes. In the network fluid regime, 
the difference between the number of bonding sites is the 




FIG. 1. Pictorial representation of a network fluid in a binary 
mixture of particles with two and tree identical patches. A 
bond between two bonding sites lowers the energy by caa- 



key parameter controlling the phase equilibria. If one 
species has at least twice the number of bonding sites of 
the other, then demixing of two network fluids gives rise 
to profound changes in the phase diagram of the system. 
The miscibility at high pressures is also controlled by this 
difference and closed miscibility gaps are present if the 
difference between the number of bonding sites is greater 
than one. 

The remainder of the paper is organised as follows. 
In section QTJ we present the model, recall Wertheim's 
theory for binary mixtures pi A[) and generalize Flory- 
Stockmayer theory of percolation to binary mixtures of 
patchy particles piBI) ; in III CI we derive expressions for 
the binary mixtures considered in this paper. In section 
lllll we present the results: phase diagrams (including per- 
colation lines) with emphasis on the network fluid regime 
and critical properties of several representative mixtures. 
Finally, in Sec. IIVI we summarize our conclusions and 
suggest lines for future research. 



II. MODEL AND THEORY 



We consider a binary mixture of N\ and N2 equisized 
hard spheres (HSs) with diameter a. The surface of each 
species is patterned with a distinct number and/or type 
of bonding sites distributed in such a manner that two 
particles can form only one single bond, involving two 
distinct sites one in each particle. A minimum distance 
between the sites is required to ensure that no sites are 
shaded by nearby bonds. A pictorial representation of 
the model under consideration can be found in Fig. [TJ 
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A. Helmholtz free energy: Wertheim's 
thermodynamic perturbation theory 

A detailed description of Wertheim's thermodynamic 
perturbation theory for pure fluids and fluid mixtures 
can be found elsewhere fill . [7J[ . Here we briefly quote the 
results and set the notation. Within perturbation theory, 
the Helmholtz free energy of a system of particles can be 
expressed as a sum of contributions from an unperturbed 
reference system where the particles interact via repulsive 
forces, here the hard-core repulsions, and a perturbation 
due to the attractive bonding interactions: 



f H = F/N = f HS + f„, 



(1) 



where N = Ni + N2 is the total number of particles, fjjs 
and ft, are the HS- and bonding- free energy per particle 
respectively, fus may be written as the sum of ideal-gas 
and excess terms: fus = fid, + fex- The ideal-gas free 
energy is given (exactly) by 



characterises the bond between a site a on a par- 



ticle of species i and a site 7 on a particle of species j. 
For simplicity, we model the interaction between sites by 
square wells with depths e ai which depend on the type 
of bonding sites (a and 7) but not on the particle species 
(i and j). As a consequence, when the species have the 

same diameter, the A^' are independent of the particle 
species, and can be written as 



_L_ 

V s 



9Hs(r) [exp(^e Q7 ) - 1] dr. (6) 



9hs ( r ) is the radial distribution function of the reference 
HS fluid and the integral is calculated over the bond vol- 
ume v a ry. We consider that all bonds have the same vol- 
ume (i.e., v ai = Vf,) and we use the contact value for the 
radial distribution function. As in previous works (20j . 
we set Vb = 0.000332285er 3 . Using these approximations 
Eq. ([6]) simplifies to 



A Q7 = — [exp(/3e Q7 ) - 1] A (r)), 



(7) 



/3/ id = lnr?-l + 



E ; 

i=l,2 



where 



(2) 



with ft = kT the inverse thermal energy, Vi the thermal 
volume and x^' = Ni/N the molar fraction of species i. 
rj — r/i + f]2 is the total packing fraction (ji = v s p, with 
p the total number density and v s = ir/6a 3 the volume 
of a single particle). The excess part, which accounts 
for the excluded volume of the monomers, is approxi- 
mated by the Mansoori-Carnahan-Starling-Leland equa- 
tion of state for HSs mixtures [l8| . The latter reduces to 
the well-known Carnahan-Starling equation of state (l9| 
when the species have the same diameter: 
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The bonding free energy is approximated by Wertheim's 
thermodynamic first-order perturbation theory [7J[. Let 
T(i) be the set of bonding sites or patches on one particle 
of species i = 1,2. The bonding free energy per particle 
is given by [Tj| 




lnXW - 



5«(T(i)) 



(4) 

where X<£ is the probability that a site of type a on 
a particle of species i is not bonded and n(F(«)) is the 
total number of bonding sites per particle of species i. 
The probabilities {Xa } are related to the total density, 
molar fractions and temperature through the law of mass 
action: 



.7=1,2 



«7 



(5) 



1 - v/2 
(1 - r?)3 



(8) 



is the contact value of gus- These approximations will 
affect the results quantitatively but are not expected 
to change them qualitatively [2(|[2l||. Substituting A Q7 
given by Eq. ([7]) into Eq. (JS{ we find that X& depends 
only on a, the type of site (i.e., Aq' = X a , V i) in 
line with the independent site approximation underlying 
Wertheim's theory. 

The equilibrium properties of the mixture are obtained 
by minimising (at a fixed composition x, pressure p and 
temperature T) the Gibbs free energy per particle g = 
p/p + fji with respect to the total density p subject to 
the constraints imposed by the law of mass action. A 
standard Newton-Raphson method is used to minimise 
g, and the law of mass action is solved simultaneously by 
a Powell hybrid method or analytically when possible. 
In what follows we will denote the composition of the 
mixture x by the molar fraction of species 1: x = x^ 
and x^ = 1 — x. 

Binodal lines are located by a standard common- 
tangent construction on g(x), which is equivalent to solv- 
ing the equations for the equality of the chemical poten- 
tials of both species in the coexisting phases (mechanical 
and thermal equilibria are satisfied by fixing the pressure 
and the temperature). 

Critical points are computed by determining the states 
which satisfy the law of mass action and the spinodal 
condition, f vv f xx — (fxv) 2 = 0. In addition, stability 
requires the vanishing of the third-order derivative in the 
direction of largest growth : 
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level i 



level i + 1 



FIG. 2. Schematic representation of two bonded particles 
(right) in a tree- like cluster (left). A particle of species j at 
level i is bonded to a particle at the previous level i—1 through 
a site a and to a particle at the next level i + 1 through a site 
P. The latter is bonded to a site 7 on a particle of species k. 



where subscripts denote partial derivatives, i.e., f xv is 
the second partial derivative of fg with respect to the 
reduced volume per particle v = I/77 and the composition 
x at constant temperature. 

The results are presented mostly through temperature- 
composition diagrams at constant pressure and pressure- 
temperature projections. The pressure, actually the os- 
motic pressure in colloidal systems, is difficult to mea- 
sure in experiments. However, temperature-composition 
diagrams at constant pressure give a general overview 
that is somewhat easier to interpret. In addition, we 
present critical temperature-packing fraction diagrams. 
From these, one obtains estimates of the packing frac- 
tions and temperature at two-phase coexistence that may 
be used to inform experimental studies. We note that the 
values of the thermodynamic variables at coexistence are 
sensitive to the model parameters (e.g., size and geom- 
etry of the molecules, volume of the bonds, interaction 
bond strength...) and thus serious attempts to predict 
the behaviour of real systems should be based on more 
accurate model building. 



level i + 1. This number is related to the number of all 
types of bonded sites on both species of particles in the 
previous level {n^} through the recursion relation (see 
Fig. [2] and reference [23[ for details) 

»8?i.7 = E E E p*-+» 

(10) 

The sum on j runs over the particle species, j = 1,2 
for binary mixtures. Td(j) is the set of different bonding 
sites on species j. f^ is the number of /3 bonding sites 
on a particle of species j and pp.^, lk is the probability of 
bonding a site /3 on a particle of species j to a site 7 on 
a particle of species k. Then, the probability of finding a 
bonded site (3 on a particle of species j is 

p ft=E E ( n ) 

k 7 er d (fc) 

which may be related to the thermodynamic variables 
through the law of mass action, since 

P Pj =l-Xf. (12) 

The equations (fTU)) are linear in {n^} anc ^ can ^ e ex ~ 
pressed in matrix form 

hi = Th Q , (13) 

where n, is a vector with components n\ k ^ and f is a 
square matrix with entries: 

T l^ = £ P Pj ^ h - Sp a ) . (14) 

The matrix T, with dimension equal the number of parti- 
cle species multiplied by the number of distinct bonding 
sites, may be diagonalized or transformed into Jordan 
form. In either case the progressions defined by Eq (fT3")) 
converge to if the largest (absolute value) of the eigen- 
values A 7fc of T is less than unity, i.e., |A 7fc < 1, 
Then, percolation occurs when \X Jk | = 1 for any value of 
Ik- 



B. Percolation in mixtures of particles with 
distinct numbers or types of bonding sites 

Recently [23j we generalised the Flory-Stockmayer 
random-bond percolation theory [l4| for a model of 
patchy particles with an arbitrary number of distinct 
bonding sites (correlated bonding probabilities). The 
theory, which neglects closed loops, was tested against 
Monte Carlo simulations and it was found to be quantita- 
tively accurate for systems of particles with three patches 
of two distinct types [24| . Here, we extend the theory and 
derive the percolation threshold for mixtures of patchy 
particles under the same no-loop assumption. 

(k) 

Consider a tree-like cluster and let n) , \ „ denote the 
number of bonded sites 7 on particles of species k at the 



C. Mixtures of particles with distinct numbers of 
identical bonding sites 

In the following we focus on the behaviour of mixtures 
of particles with distinct numbers of identical bonding 
sites, say A. There is then a single bonding energy €aa, 
which sets the energy scale. We will denote by f]/ — f A 2 ^ 
a binary mixture of particles of species 1 characterised by 
/j^ (A) sites and particles of species 2, characterised by 

/f (A) sites. _ 

For these mixtures the bonding free energy given by 
Eq. ((3]) simplifies to 

Pf b ={f)Ux A -^+±y (15) 
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where 



kT/8 



(f)=xfV + (l-x)f% ! , 



(2) 



(16) 



is the average functionality or average number of bonding 
sites per particle at a given composition x. The bonding 
free energy can be split into two terms: the bonding en- 
ergy Ub and an entropic term related to the number of 
ways of bonding two particles. As the sites are indepen- 
dent, Ub is simply, 



/3 



Ub 
N 



1-X A 



(f)PtAA, 



(17) 



while the entropic term is obtained by subtracting the 
total bonding free energy, given by Wertheim's theory, 
from Ub- 

The law of mass action given by Eqs. ([5]) reduces to a 
single equation for the fraction of unbonded sites, 

1 - X A = x/^vX^Aaa + (1 - x)ffr,X\ A AA , (18) 

and can be solved analytically. The percolation matrix 
is a 2 x 2 square matrix with entries 

if — ( PA^A, (fA } - 1) PA 2 ^A, (f A 2) - 1)| (1Q) 

' \PA^A 2 (f A 1} 1) PA^AMT - iy ' 

The probabilities p Ai ^ A . are found by relating the prob- 
ability of finding a bonded site A, 

Pa = PA X ^A X + PA X ->A 2 = PA 2 ^Ai + PA 2 ^A 2 , (20) 

to the fraction of unbonded sites P A = 1 — X A . A term- 
by-term analysis of Eqs. (|18|) and (|2U|) gives 

p Al ^ Al =p A2 ^ Al = xfx\x A A AA , 
PA 1 - f A 2 = Pa 2 ->a 2 = (1 - x)f A 2) i]X A A AA . (21) 
The only non-zero eigenvalue of T is 

A = p A2 ^ A2 (f A 2) - 1) +PA 1 ->A 1 (f { A ) - 1), (22) 

and thus the system is percolated if A > 1 . The probabil- 
ity of finding a bonded site at the percolation threshold, 
P p , can be found by setting A = 1: 



P, 



xf^ + (1 - x)f\ 



(23) 



Note that this probability depends only on the function- 
ality of the particles and on the composition of the mix- 
ture. Furthermore it reduces to the Flory-Stockmayer 
result in the limit of pure fluids: xW = 1 P p = 

l/Uf ~ !)• 



III. RESULTS 

We begin by discussing briefly the properties of the 
pure fluids as they will be used in the following sections 
in the analysis of the phase behaviour of the mixtures. 
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FIG. 3. Pressure (a), packing fraction (b) and fraction of un- 
bonded sites (c) as a function of temperature in pure fluids of 
particles with 3, 4 and 5 identical bonding sites. Circles indi- 
cate critical points. The dotted horizontal line in (a) marks 
the pressure p* = pv s /eAA = 1.047 x 10~ 4 referred to in the 
text. Percolation lines are dashed in panels (a) and (b). The 
system is percolated above the percolation lines. In panel (c) 
the percolation thresholds (omitted for clarity) are horizontal 
lines with values Xa = 1 — — 1), being /a the number 

of patches. 



A. Pure fluids 

In Fig [3] we summarise the properties of pure fluids of 
particles with 3, 4 or 5 identical bonding sites. These 
fluids were previously studied using Wertheim's theory 
and Monte Carlo simulations [ll|. At temperatures be- 
low the critical temperature the system undergoes a first- 
order phase transition as pressure increases. The tran- 
sition involves two fluid phases with different densities 
and fractions of unbonded sites. In what follows we de- 
note these phases by liquid L (higher density and smaller 
fraction of unbonded sites) , and vapour V (lower density 
and larger fraction of unbonded sites). The LV coexis- 
tence line ends at a critical point which moves towards 
lower pressure and temperature as the number of bond- 
ing sites decreases. The packing fraction and the frac- 
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tion of unbonded sites at the critical point also decrease 
as the number of bonding sites decreases. As the tem- 
perature vanishes, the packing fraction of the coexisting 
liquid phase saturates to a value that decreases as the 
functionality decreases. The stability of the liquid phase 
increases as the number of bonding sites per particle in- 
creases but otherwise the phase behaviour is qualitatively 
the same for all systems. 

The percolation line (dashed lines in Fig. [3]) intersects 
the coexistence line below the critical point on the vapour 
side. Therefore, the liquid phase is always a network 
fluid. Near the critical point two percolated states or 
network fluids coexist, confirming that percolation is a 
pre-requisite for criticality in systems where the bonding 
interactions are attractive [25| . 

The behaviour is completely different for a system of 
particles with 1 or 2 bonding sites. In either case there 
is no LV phase transition. Systems with 1 or 2 bond- 
ing sites will form dimers and linear chains, respectively. 
The absence of clustering of the dimers and of branch- 
ing of the chains prevents the fluids from condensing. A 
polymerization transition occurs for linear chains in the 
limit of vanishing temperature, T — > [26[. 



B. Binary mixtures of particles with distinct 
functionality: 2a — 3a mixture 

We start by considering a 2a — 3 a mixture (i.e., a bi- 
nary mixture where species 1 has two patches of type 
A and species 2 has three). The critical properties of 
this system were investigated recently by grand-canonical 
Monte Carlo (MC) simulations [lfj- The results indicate 
that the critical packing fraction and temperature de- 
crease continuously towards zero as the fraction of 2a 
particles in the mixture approaches one (i.e., the average 
functionality decreases towards two). The authors of jlOj ] 
mapped the properties of the binary mixture to those of 
a pure fluid of particles with a (non-integer) number of 
bonding sites equal to the average functionality (/). Al- 
though the critical parameters were predicted correctly 
by this mapping, the effect of the entropy of mixing on 
the phase behaviour of the mixture was overlooked, and 
this, together with a full description of the phase dia- 
gram, will be addressed in this section. 

In Fig. 2]we illustrate the results for the phase diagram 
of the 2^ — 3,4 mixture in the T— x plane at constant pres- 
sure. The diagrams correspond to three distinct pressures 

(3) 

below p c , the critical pressure of the pure 3,4 system. 

(3) 

The mixture is always stable for pressures p > p c (we 
will return to this point later). At lower pressures the 
LV phase transition of the pure 3a system shifts from 
the x = axis to finite values of the composition, as 2a 
particles are added to the mixture. The two-phase LV 
region ends in a lower critical point at (x c ,T c ) and in- 
creases as the pressure decreases. The percolation line 
intersects the binodal on the vapour side, close to the 
critical point. Below the percolation line the system is a 
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FIG. 4. Phase diagrams (at different pressures) in the scaled- 
temperature vs composition (x = x^ 1 ') plane of a binary 
mixture of particles with 2 (species 1) and 3 (species 2) 
identical bonding sites, (a) pv s /eAA = 5.236 x 10~ 4 . (b) 
p* = pv s /eaa = 1.047 x 1CT 4 . (c) pv s /eaa = 1.885 x 1(T 7 . 
The shaded area indicates the two-phase region. Circles de- 
note critical points. Percolation lines are dashed. The insets 
depict scaled-temperature vs packing fraction plots along the 
binodal curve. 



network fluid in the sense that there is a non-zero proba- 
bility of finding a (more or less transient) infinite cluster. 
Coexistence involves two network fluids in a finite range 
of temperatures above the critical temperature, which 
decreases as the pressure decreases. 

An increase in the fraction of 2a particles reduces 
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FIG. 5. Cut of the (T, rj, x) phase diagram at constant compo- 
sition for a binary mixture of particles with 2 (species 1) and 
3 identical bonding sites. Shaded areas indicate two-phase 
regions. Open circles denote critical points. Percolation lines 
are dashed: Black and brown for mixtures with x = 0.60 and 
x — 0.95, respectively. 
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FIG. 6. Critical properties of a pure fluid with non-integer 
number of bonding sites (red dotted lines) and a binary mix- 
ture of particles with 2 and 3 identical bonding sites (black 
dashed lines). Grey circles mark the points corresponding to 
a fluid of particles with 3 identical sites, (left) Scaled-critical 
temperature T* = kT/eAA vs critical packing fraction. Mid- 
dle and right panels depict the scaled-critical temperature and 
critical packing fraction as a function of (/) , the average num- 
ber of bonding sites per particle at the critical point (see text). 



the probability of branching, which drives condensation. 
Therefore, a decrease in the stability of the liquid phase 
is expected as x increases, and, at a given pressure, the 
vapour is stable at temperatures below the transition 
temperature of the pure 3,4 system (the temperature at 
the intersection of the binodal with the x = axis). How- 
ever, inspection of the insets of Fig. @] (T — rj plots along 
the binodal curves) reveals that the density of the liquid 
phase decreases monotonically with T down to T c . This 
means that the reduced stability of the liquid phase in the 
T — x plane is accompanied by an increase in the range 
of densities where the liquid phase is stable. In order to 
show this, a cut of the (T, rj, x) phase diagram at constant 
composition is depicted in Fig. [5] It is important to note 
that, in this representation, a tie line connecting two co- 
existing points in the binodal is out of the plane, except 
for the critical points and the pure fluid at x = 0. The 
two-phase region (shaded area) extends over the whole 
range of composition (except x = 1), and decreases as 
the composition of 2a particles increases. At x = 0.95, 
for example, it is possible to stabilise a liquid phase with 
rj < 0.009 while the lowest density of the liquid in a pure 

3a system is rj^ ~ 0.09 (see Fig. [3J (b)). In fact, as 
shown by Sciortino et al. [10j . r\ c and T c approach zero 
asymptotically as the pressure vanishes. 

The critical line of the mixture is depicted in Fig. [6] 
where it is compared to the line of critical points of a 
pure fluid with a number of bonding sites equal to the 
average functionality of the mixture at the critical point 
[Tol | . As expected, there are differences between the criti- 
cal properties of the mixture and those of the pure fluids. 
The differences result from the composition fluctuations 
in the mixture, neglected in the mapping to a pure fluid 



[27|- Nevertheless, the mapping describes correctly the 
most salient feature of the 2 a — 3a mixture critical line: 
The vanishing of the critical parameters as the composi- 
tion of the 2a — 3a approaches x = 1 . 

Indeed, the agreement between the theoretical predic- 
tions for the critical point of pure fluids with average 
functionality and the MC simulations of the critical line 
of 2a — 3a mixtures is quite remarkable (see Fig. 3. of 
reference [loj]). Nevertheless, the structure of the criti- 
cal line revealed by the MC simulations signals the effect 
of composition fluctuations, confirmed by the theoretical 
calculations of the mixture's critical line (left panel of 
Fig. [5]). In addition, both the MC and the theoretical 
results for the critical parameters T c and rj c of the bi- 
nary mixture are above those of the pure fluid with the 
same functionality, as illustrated in the middle and right 
panels of Fig. [5] 



C. Binary mixtures of particles with distinct 
functionalities 

In this section we consider the binary mixtures: 2a — 
4a, 2a — 5a and 1a — 3a- The results are summarised in 
the next three figures. Temperature-composition phase 
diagrams at constant pressure are depicted in Fig. [7] 
Results are shown at different pressures below and above 

(3) 

Pc . In all cases the phase diagram at pressure p\ is de- 
picted (the phase diagram of a 2a — 3a binary mixture at 
the same pressure is plotted in panel (b) of Fig. |3|) . The 
properties of the critical lines are analysed in Fig. [SI and 
pressure-temperature projections of the phase diagrams 
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are shown in Fig. 



1. 2a — 4a mixture: liquid- liquid demixing driven by the 
entropy of bonding 

In panel (a) of Fig. [JJwe plot the T—x phase diagram of 
2a — 4a mixtures at different pressures below and above 
pif^ (the critical pressure of the 4 a fluid). At intermedi- 
ate pressures (for example p*) the T — x phase diagrams 
are qualitatively the same as those of the 2a — 3 a mix- 
ture: A two-phase region bounded by a lower critical 
point starts at the x = axis and the percolation line in- 
tersects the binodal on the vapour side, near the critical 
point. However, by contrast to the 2a — 3 a mixture, the 
LV density gap along the binodal curve (not shown) has 
a non-monotonic behaviour. It increases near the x = 
axis and then decreases until it vanishes at the critical 
point. This implies that the drive for phase separation 
increases when a small fraction of 2a particles is added to 
the pure 4a fluid, which has not been observed in the 3a 
case. This behaviour may be understood in terms of the 
balance between the entropy of mixing and the entropy of 
bonding. Consider a bond between particles with 2 and 
n sites: The resulting structure is a two-particle cluster 
with n sites available for bonding. The structure that 
results from a bond between two particles with n sites 
is a two-particle cluster with 2(n — 1) sites available for 
bonding. The latter has n — 2 additional sites that are 
available for bonding, and thus the loss in the entropy 
of bonding increases as n increases, while the gain in the 
entropy of mixing remains the same. This suggests that 
the tendency for phase separation increases with n, as 
observed in the 2a — 4a mixture. 

A more striking consequence of a stronger drive for 
phase separation is the behaviour of the 2a — 4a mixture 
above the critical pressure of the pure 4a fluid. As shown 
in panel (b) of Fig. [9] there is a range of pressures above 
p^ where phase separation still occurs. The two-phase 
region is a closed loop bounded above (below) by an up- 
per (lower) critical point. An example of the T — x phase 
diagram in this range of pressures is plotted in panel (a) 
of Fig. [7J (darkest grey). The percolation line is always 
above the phase separation region. As the pressure in- 
creases the demixing region decreases, and it vanishes at 
a given pressure. In panels (a) and (b) of Fig. [H] we com- 
pare the pressure-temperature projections of the phase 
diagrams of 2a — 3 a and 2a — 4a mixtures. The critical 
line has a maximum at p > pi 4 ^ for the 2a — 4a mix- 
ture which is absent in the 2a — 3a case. As a result, 
closed miscibility gaps are found in a range of pressures 
delimited by this maximum and p^ ■ 

Based on the topology of the pT projections, both the 
2a — 3a and 2a — 4a mixtures are limiting cases of type I 
mixtures (accor ding to the classification of van Konynen- 
burg and Scott J28|) where one of the species (particles 
with 2 bonding sites) has vanishing critical temperature. 



In type I mixtures the critical line is continuous and it 
connects the critical points of the pure fluids. In 2a — 3a 
and 2a — 4a mixtures the critical line starts at the criti- 
cal point of species 2 and tends asymptotically to T = 
as the pressure vanishes, p —> (species 1 undergoes a 
polymerization transition as T — > 0) . From a topological 
point of view the significant difference between the crit- 
ical behaviour of 2a — 3a and 2a — 4a mixtures is that 
the critical line is monotonic in the former while it is non- 
monotonic in the latter, with closed miscibility gaps near 
the critical region of the 4a fluid. However, as we will see 
now, there is another important difference between both 
mixtures at very low pressures. 

Let us focus on the temperature-composition phase di- 
agram of a 2a — 4a mixture at very low pressure (panel 
(a) of Fig. [JJ medium gray) . We find a large two-phase 
region, the shape of which suggests the presence (or prox- 
imity) of two different phase transitions. Near x = 0, 
we observe a LV transition also present in the 2a — 3a 
mixture, but as the temperature is lowered, a new type 
of demixing appears (as a bulge) in the phase diagram. 
In this region two percolated or network fluids coexist 
(note that this region is well inside the percolated area). 
We will refer to this as liquid-liquid coexistence in order 
to distinguish it from the LV coexistence. Both regions 
form a continuous two-phase region but their different 
origin is revealed by analysing the critical points (see Fig. 
[5]). As the pressure decreases, the critical temperature 
vanishes but the critical density tends to a finite value, 
r/ c — > 0.053. The average functionality at the critical 
point approaches asymptotically (/) = 2.18 (x c = 0.911) 
rather than (/} = 2 (x c = 1). Therefore, although the 
topology of the mixture is unaffected, the critical line 
changes its character in a continuous fashion from liquid- 
vapour (near the critical point of the 4 patches fluid) to 
liquid-liquid. The new liquid-liquid transition preempts 
the liquid- vapour phase transition and there is no liquid- 
liquid-vapour triple line in this mixture. Nevertheless, 
LLV triple lines are observed in other mixtures as we 
will discuss later. 

In Fig. [8] we have also presented the critical properties 
of pure fluids with a number of bonding sites equal to 
the ave rag e functionality of the mixtures at the critical 
points [lGj . It is clear that a mapping of the properties 
of the mixture to those of a pure fluid with the same 
functionality completely fails in this case, as no empty 
liquids are found at low pressure. 

Our results strongly suggest that the new liquid-liquid 
demixing results from a balance between the entropy of 
bonding and the entropy of mixing. Let us consider a 

(2) 

bond between two particles with f A sites. The result- 

(1) 

ing structure is a two-particle cluster with 2{f A ' — 1) 
unbonded sites. This corresponds to an increase of 

Ia^ ~ fA unbonded sites when compared to the struc- 
ture that results from bonding particles of species 1 and 
2. When this increase is large compared to the num- 
ber of bonding sites of species 1, the entropy of bonding 
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FIG. 7. Phase diagrams at different pressures for mixtures of particles with distinct numbers of identical bonding sites in the 
scaled-temperature kT/tAA vs composition x plane. Shaded areas are two-phase regions, open circles denote critical points and 
percolation lines are dashed, (a) 2 a — 4a mixture: particles with two A-sites (species 1) and particles with four A-sites (species 
2), pressures: pv s /eAA = 2.513 x 10~ 3 (darkest grey), pl — pv s /tAA = 1.047 x 10 -4 (lightest grey) and pVs/tAA = 1.885 x 10~ 7 
(medium grey), (b) 2a — 5a mixture, pressures: pVs/^AA = 1.047 x 10 -2 (dark grey), p — pl (light grey), (c) 1a — 3a mixture, 
pressures: pv a /^AA = 1.047 x 10 -3 (dark grey), p = p\ (light grey), (d) 3a — 4a mixture, pressures: pVs/^AA = 1.047 x 10~ 3 
(dark grey), p = pl (light grey) . The inset in (d) is a zoom of the region close to the LV critical point at pl. 



gained by bonding particles of species 2 compensates the 
entropy of mixing lost and drives the liquid-liquid phase 
separation. Considering only an integer number of bond- 
ing sites the condition can be written as: 

/f > 2/«. (24) 

That is, if one species has at least twice the number of 
bonding sites of the other, then liquid-liquid demixing, 
driven by the entropy of bonding, occurs. Nevertheless, 
if we allow the particles to have a non-integer number 



of bonding sites, then the condition for finding liquid- 
liquid demixing is not exactly given by Eq. (|2"4"|) . For 
example, in mixtures where species 1 has two sites, we 
have found empty liquids - and no liquid-liquid demixing 
- if < 3.99, while liquid-liquid demixing occurs when 

f { 2 ] ^ 3.99. m. 

We have checked the heuristic argument of Eq. (pH)) 
by calculating the phase diagrams of other mixtures sat- 
isfying this condition. 
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FIG. 8. Critical properties of binary mixtures of particles 
with distinct numbers of identical bonding sites, (a) Scaled- 
critical temperature kT c /eAA vs critical packing fraction r) c . 
(b) Scaled-critical temperature as a function of the average 
number of bonding sites at the critical point (/). The be- 
haviour of a hypothetical pure fluid with non-integer number 
of sites is also plotted as a red-dotted line. Grey circles denote 
critical points of the pure fluids with 3, 4 and 5 sites. 
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FIG. 9. Pressure-temperature projections of the phase dia- 
grams of different binary mixtures, (a) 2 a — 3a mixture, (b) 
2a — 4:A mixture, (c) 2a — 5a mixture, (d) 3a — 4a mixture. 
Solid curves are the liquid-vapour transition lines of the pure 
fluids. Dashed lines are the critical lines of the mixture. Open 
circles denote the critical points of the pure fluids. 

(/) — > 3.22). There is a slight difference in the behaviour 
of 2a — 4:A and 2a — 5a mixtures. In the latter, the 
critical packing fraction is always higher than the crit- 
ical packing fraction of the pure fluid with 5 patches, 

(5) 

rjc , and increases monotonically with decreasing pres- 
sure (see panel (a) of Fig. [5]). In the former (2a — ^a 
mixture) , the critical packing fraction increases from r/i 3 ^ , 
reaches a maximum, then decreases slowly, and tends to 
a plateau that is lower than the critical density of the 
pure fluid with 3 patches. 



2. 2a — 5a mixture 

Phase diagrams of the 2a — 5a mixture are depicted in 
panel (b) of Fig. [7] at two different pressures below (light 

grey) and above (dark grey) Pc 5 \ the critical pressure 
of the pure fluid with 5 bonding sites. As in the pre- 
vious mixture, there is a region of pressures above pi 
where closed miscibility gaps are present. In this high 
pressure region, the phase coexistence involves two per- 
colated fluids. As the pressure decreases, the two-phase 
region grows very rapidly, occupying a large fraction of 
the phase diagram (compare, for example, the phase dia- 
gram of this mixture at pressure pi with other mixtures 
at the same pressure). 

Topologically this mixture is still a limiting case of type 
I. The analysis of the percolation threshold and the crit- 
ical properties (see Fig. [5]) indicates the coexistence of 
two network fluids in the limit of low pressures. As the 
pressure decreases, the critical temperature tends asymp- 
totically to T c = 0, but the critical packing fraction and 
critical composition tend to finite values (tj c — >• 0.282 and 



3. 1a — 3a mixture 

We proceed by considering the 1^ — 3,4 mixture (T — x 
phase diagrams are plotted in panel (c) of Fig. [7]) which 
also satisfies Eq. (|2~4"| . The addition of 1a particles to 
the 3a fluid drastically reduces both the probability of 
branching and the average cluster size. At x as 0.75 the 
predominant structures are isolated clusters of three 1a 
particles that saturate the bonds of the 3 a particles. The 
most obvious consequence of this is that when T — > 0, 
the percolation line tends asymptotically to x = 0.75 
rather than to x = 1. Thus, it is not possible to observe 
a percolated fluid if the composition of the mixture is 
greater than x = 0.75. 

We also found closed miscibility gaps at p > p c \ 
an example is the T — x phase diagram at pv s /eAA = 
1.047 x 10~ 3 depicted in panel (c) of Fig. [7J(dark grey). 
The line of critical points, represented in Fig. behaves 
qualitatively as that of the 2a — 5a mixture and, sur- 
prisingly, it is possible to find coexistence between fluid 
phases with an average functionality smaller than two 
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(note that (/) -> 1.98 when T c -> 0). Topological^, this 
is still a type I mixture. 

(2) 

D. From type I to type V mixtures: 3a — f\ 
mixtures 

(2) 

Finally, we focus on 3a — f A mixtures, in particular, 
3 a — 4a and 3 a — 6a- The latter satisfies the liquid- liquid 
demixing condition, Eq. ([HJ) , while the former does not. 
By contrast to the mixtures considered previously, both 
pure fluids exhibit LV transitions. 

The phase diagram of the 3,4 — 4 a mixture is plotted 
in panel (d) of Fig. [7] At intermediate pressures, be- 
tween p c 3 ^ and p c ^ , the phase diagram consists of a small 
two-phase region bounded by a lower critical point. The 
percolation line starts at x = 1 at finite temperature and 
intersects the binodal on the vapour side, close to the 
critical point. At pressures below p c (e.g., p\) the criti- 
cal point disappears and there is LV coexistence over the 
whole range of composition. The liquid phase is always 
percolated in this range of pressures. Above p c ^ the sys- 
tem is completely miscible. The line of critical points in a 
pT projection, represented in panel (c) of Fig. |H1 connects 
the critical points of the two pure fluids and thus it is a 
standard type I mixture. The temperature and the den- 
sity on the critical line are similar to those obtained from 
the mapping to a pure fluid with the same functionality 
(see Fig. [8]). In 3a — 5a mixtures, the only difference 
is the existence of closed miscibility gaps at pressures 
higher than p c ^ (the critical line is non-monotonic. See, 
for example, the pT projection of this mixture in panel 
(d) of Fig. 

However, new phenomenology is found in the phase 
behaviour of 3a — 6a mixtures. Three representative 
phase diagrams are plotted in Fig. [TO] At pressures 
slightly below p c 6 ' there is a large two-phase region start- 
ing at x = and ending in a lower critical point at fi- 
nite values of the composition (panel (a) of Fig. HU|) . 

(3) 

As p c is approached, a phase transition involving two 
liquid-like states is clearly visible. There are two criti- 
cal points, where liquid-liquid and liquid-vapour coexis- 
tence end (see panel (b)). A liquid- liquid- vapour triple 
point is present at temperatures slightly above the tem- 
perature of the liquid-vapour critical point. Below p c 
the liquid- vapour critical point disappears but the liquid- 
liquid critical point is present down to p^ w 0.43pc 3 ^- 
At pressures below p(' the liquid-liquid demixing be- 
comes metastable with respect to the liquid- vapour tran- 
sition (see panel (c)) and the phase diagram is similar to 
those of the 3a — 4a or 3a — 5a mixtures. 

The critical properties of mixtures where species 1 has 
three bonding sites are shown in Fig. [TT](a). For 3a — 4a 
and 3a — 5a mixtures the critical line is continuous and 
connects the critical points of both fluids. In the 3a — 6a 
mixture, however, the critical line is discontinuous (see 
also the pT projection in panel (b) of the same figure). A 
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FIG. 10. Phase diagrams (at different pressures) in the scaled- 
temperature vs composition (x = a;' 1 ') plane of a binary mix- 
ture of particles with 3 (species 1) and 6 (species 2) identical 
bonding sites (3a — 6a mixture), (a) pv s /eAA = 5.236 ■ 10~ 3 . 
(b) pv s /eAA = 1.047 • 10" 3 . (c) pi = pv s /€ A A = 1.047 • HP 4 . 
The shaded area is the two-phase region. Circles denote the 
critical points. Percolation lines are dashed. The horizontal 
dotted line in (b) indicates the position of the triple point. 
The inset in (b) is a zoom of the region near the x = 1 axis. 
In this zoom the percolation line connects the x = 1 axis 
to the binodal. The dashed-dotted line inside the demixing 
region in (c) is the metastable binodal for LL coexistence. 
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FIG. 11. (a) Scaled-temperature vs packing fraction at the 
critical lines of various 3a — f\ mixtures. The critical line of a 
hypothetical pure fluid with a non-integer number of bonding 
sites is also represented as a red-dotted line. Grey circles 
denote the critical points of pure fluids with 3, 4, 5, and 6 
bonding sites. Blue-squares are the critical end points of the 
3a — 6a mixture, (b) Pressure-temperature projection of the 
phase diagram of the 3a — 6a mixture. Solid curves ending in 
an open circle are the LV lines of the pure fluids. The other 
solid line (see zoom in the inset) is the LLV three phase line. 
The critical lines are dashed. 



liquid-vapour critical line extends from the critical point 
of the component with a lower number of bonding sites 
(3a) and terminates at a three phase line at an upper crit- 
ical end point. Another liquid- vapour critical line starts 
at the critical point of the other component (6^), changes 
its character continuously to a liquid-liquid critical line 
and ends on the same three phase line at a lower critical 
end point (type V binary mixture). 

The same topological change occurs in mixtures where 
species 1 has 4 bonding sites (see Fig. [12"]) when species 
2 has 8 bonding sites. 

We have not studied mixtures with > 4, as the 
predictions of Wertheim's theory become less accurate 
as the number of bonding sites increases [3(j, and also 
because the density is not sufficiently low to rule out the 
presence of stable solid phases that will appear at high 
pressure and/or low temperature. However, the stabil- 
ity of the solid is very sensitive to the location of the 
bonding sites: asymmetric distribution (e.g., tetrahedral 
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FIG. 12. Scaled-temperature vs packing fraction at the crit- 

(2) 

ical points of 4a — f\ mixtures. The behaviour of a hypo- 
thetical pure fluid with a non-integer number of bonding sites 
is also represented as a red-dotted line. Grey circles denote 
the critical points of pure fluids with an integer number of 
bonding sites in the range 4 to 8. Critical end points of the 
4a — 8a mixture are denoted by blue-squares. 



symmetry in particles with 4 bonding sites) will promote 
the formation of crystalline phases, whereas the oppo- 
site occurs if the bonding sites are randomly distributed. 
In Wertheim's theory the bonding sites are independent, 
and thus it will describe more accurately particles with 
bonding sites that are distributed randomly. 



IV. DISCUSSION AND CONCLUDING 
REMARKS 

We have analysed the phase behaviour, percolation 
threshold and critical properties of binary mixtures of 
patchy colloidal particles, using Wertheim's theory and 
an extension of the Flory-Stockmayer theory of percola- 
tion for binary mixtures. We have restricted the study 
to mixtures where all bonding sites are identical (there is 
a single bonding probability) and the particles have the 
same diameter (no drive for phase separation from the 
excluded volume of the monomers) . Despite the simplic- 
ity of the model, we have found a rich phase behaviour 
including closed miscibility gaps, liquid-liquid phase sep- 
aration and topological changes in the phase diagram, as 
the number of bonding sites of the particles is varied. 

We have found that the entropy of bonding plays a cru- 
cial role in the stability of fluid phases and drives novel 
liquid-liquid phase transitions, not present in mixtures 
of simple fluids. Consistent with this observation the 
closed miscibility gaps and the liquid-liquid demixing oc- 
cur deep in the percolated region of the phase diagram. 
It is then two network or structured fluids that coexist 
at these new phase transitions. 

The difference in the number of bonding sites of both 
species was found to be the key parameter that controls 
the topology of the phase diagram and the nature and 
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number of the fluid phase transitions of the mixture. In 
particular, we found that: 

1. The phase diagram is type I (the critical line of 
the mixture connects the critical points of the pure flu- 
ids) if the species with a higher number of bonding sites 
has less than twice the the number of bonding sites of 
the other species. If this condition is not satisfied, then 
the entropy of bonding drives a new liquid-liquid phase 
separation and the mapping to a pure fluid with average 
functionality fails. Furthermore, if both species have fi- 
nite LV critical points the phase diagram changes from 
type I to type V (A LV critical line extends from the criti- 
cal point of the fluid with the smaller number of bonding 
sites and terminates at a three phase line at an upper 
critical end point. Another LV critical line starts at the 
critical point of the other species, changes its character 
continuously to a liquid-liquid critical line and ends on 
the same three phase line at a lower critical end point). 

2. The critical line is monotonic if the difference be- 
tween the number of bonding sites of the two species is 
one: 

\f ( A ) -fA\ = l- (25) 

This is the case for 2^—3,4, 3,4—4,4 and 4^ — 5^4 mixtures, 
and the corresponding fluids are completely miscible at 
all pressures. The mapping to a pure fluid with average 
functionality is (qualitatively) correct. By contrast, the 
critical line is non-monotonic if the difference between 
the number of bonding sites of the two species is greater 
than one: 

\f A 1] -f A 2) \>l. (26) 

Examples are 1^ — 3^, 2^ — 4^, 3a — 5 a or 4a — 6a mix- 
tures. In this case closed miscibility gaps appear at high 
pressures and the mapping to a pure fluid with average 
functionality fails. 

Changes in the topology of the phase diagram of bi- 
nary mixtures have been reported experimentally and 
theoretically for a variety of systems. For example, the 
CO2 + n— alkane homologous series exhibits transitions 
from type I to type II, IV and III as the chain length of 
the hydrocarbon increases [3l|. These changes are gen- 
erally correlated to the ratio of critical temperatures and 
pressures of the pure fluids. In the binary mixtures dis- 
cussed here these ratios also vary when the number of 
bonding sites varies. We believe, however, that the driv- 



ing force for the change in this class of mixtures is the dif- 
ference in the bonding entropy associated with bonds be- 
tween like- and unlike- particles. It is possible to address 
this question by considering mixtures of particles with 
different diameters (or different types of bonding sites) 
to constrain both fluids to have similar critical pressures 
and temperatures. 

Related studies on the global phase diagrams of binary 
mixtures using Wertheim-based equations of state can be 
found in the literature. See, for example, Ref. jjj^, where 
binary mixtures of water and n— alkanols are modeled 
by the statistical associating fluid theory (for a recent 
review see, for example, Ref. [33]). In these models, 
however, the attractive interactions are not limited to 
the interaction between bonding sites as considered in 
the present study. This limit is unlikely to be relevant 
for molecular fluids but appears to be relevant to the 
newly synthesized patchy particles, and gives rise to new 
and interesting phenomenology. 

We have already studied the empty fluid regime of bi- 
nary mixtures of particles with 2 and 3 patches of differ- 
ent types [U with the goal of designing structured fluids, 
with novel macroscopic properties, including stable bigel 
phases (35| . 

As a final remark, we note that a state-dependent func- 
tionality / may be defined by mapping the moments 
of the cluster distributions functions of the mixture, to 
those of pure fluids, similarly to what is done in (24|. A 
pure fluid characterised by a state dependent functional- 
ity may be useful to describe aspects of the mixture phase 
behaviour, such as closed miscibility gaps or liquid-liquid 
phase separation, that are not described by the mapping 
to a fluid with average functionality. 
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